Near-threshold high-order harmonic spectroscopy with aligned molecules 
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We study high-order harmonic generation in aligned molecules close to the ionization thresh- 
old. Two distinct contributions to the harmonic signal are observed, which show very different 
responses to molecular alignment and ellipticity of the driving field. We perform a classical elec- 
tron trajectory analysis, taking into account the significant influence of the Coulomb potential on 
the strong-field-driven electron dynamics. The two contributions are related to primary ionization 
and excitation processes, offering a deeper understanding of the origin of high harmonics near the 
ionization threshold. This work shows that high harmonic spectroscopy can be extended to the 
near-threshold spectral range, which is in general spectroscopically rich. 

PACS numbers: 



High-harmonic generation (HHG) in strong-field light- 
matter interactions is generally understood in terms of a 
three-step process: an electron is removed from the tar- 
get, accelerated, and driven back to its parent ion by 
the laser field. Finally, it recombines into the bound 
state, and emits a high-harmonic photon [l|. This pro- 
cess has been extensively studied for harmonics with en- 
ergies higher than the ionization threshold (ip), corre- 
sponding to frequencies in the extreme ultraviolet (XUV) 
range. In this letter, we focus on the intricate region close 
to the threshold and use aligned molecules to elicit the 
dynamics of near-threshold HHG. 

The spectral region near the ionization threshold is, in 
general, spectroscopically very rich. In this respect, using 
near-threshold harmonic spectroscopy to probe specific 
structures such as resonances is an appealing alternative 
to standard photoionization spectroscopy. It could help 
uncover nontrivial attosecond dynamics occuring near 
the threshold, such as excitation of electronic wavepack- 
ets or ionization via resonant states 2, 3]. However, the 
extension of high harmonic spectroscopy to the thresh- 
old region requires a detailed knowledge of the underlying 
HHG mechanisms in this specific regime. 

Furthermore, high harmonic spectroscopy is being used 
in a growing number of experiments to measure the struc- 
ture and dynamics of small atomic and molecular systems 
with unprecedented sub-femtosecond and Angstrom res- 
olution l^-lSl- As we move to more complex molecular 
structures, several orbitals contribute to the laser-target 
interaction and nontrivial hole dynamics are induced. 
Due to the contribution of multiple orbitals, there is no 
clear distinction between near- and above-threshold har- 
monics. The ability to identify the different processes in 
this regime is an essential step towards decoupling and 
identifying the contributions from the various orbitals. 

Above the ionization threshold, the intuitive three-step 
model is qualitatively confirmed and quantitatively rein- 
forced by the Strong Field Approximation (SFA)[10(. In 
this model, a photon at frequency u is emitted in the 



energy range I p < hui < I v + 3.17 • U p , where I p is the 
ionization potential and U p the ponderomotive energy. 
In the SFA, the free electron is controlled by the strong 
laser, whereas the Coulomb potential is neglected. The 
SFA predicts many of the fundamental characteristics of 
HHG, in particular, the fact that each harmonic mainly 
stems from two electronic pathways, commonly referred 
to as "short" and "long" trajectories [llj]. However, the 
SFA is not accurate as we approach the ionization thresh- 
old, where the Coulomb force cannot be neglected [12 ) . 

There has recently been increased interest in the near- 
threshold regime [l3Ml5j . Yost et al. have considered 
the generation of harmonics below the ionization thresh- 
old in xenon atoms [l3j . In conjunction with measure- 
ments of harmonic yields as functions of the laser inten- 
sity, they solved the time-dependent Schrodinger equa- 
tion and identified two distinct contributions to below- 
threshold harmonics, with different intensity-dependent 
phases. The contribution with the highest intensity de- 
pendence of the phase was attributed to long electron tra- 
jectories, whereas the second contribution, with smaller 
intensity dependence, was related to primary multipho- 
ton processes. Hostetter et al. [HI supported this obser- 
vation from a semiclassical point of view. 

How can we directly and experimentally distinguish 
between the different HHG processes in the vicinity of 
the ionization threshold? There are two important knobs 
that can be used to manipulate the HHG process. The 
first knob is related to the structure of the bound state 
and can be manipulated by producing harmonics from 
aligned molecules. The second is related to the free elec- 
tron dynamics and can be tuned by controlling the laser 
ellipticity. Here we combine the two knobs, using exten- 
sively studied N2 and O2 molecules, and manipulate both 
the target structure, through alignment, and the laser 
field, which guides the free electron wavepacket. Such a 
manipulation is applied to elicit the dynamical features 
of threshold HHG. In addition, we perform a theoretical 
analysis in terms of classical electron trajectories, which 
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FIG. 1: Harmonic yields in O2 aligned at (left) and 90 
(right) degrees, for long (bottom) and short (top) trajectories, 
as functions of the ellipticity e of the probe pulse. The signal 
from each harmonic is normalized to its value at e = 0. The 
inset shows the full scale of the increase in harmonic 9. 



yields an intuitive picture of the interaction and improves 
our understanding of the underlying mechanisms. 

The experiment is composed of two parts. We con- 
trol the molecules using a pump-probe setup in which 
the pump creates rotational w avep ackets and the probe 
generates the harmonic signal [16j . Once the molecular 
axis is fixed we manipulate the free electron's dynam- 
ics by scanning the ellipticity of the probe beam. High 
harmonics are generated with an 800-nm, 30-fs, 1-KHz 
laser beam, focused into a pulsed gas jet at ~ 1.5 x 10 14 
W/cm 2 intensity. As a pump we use part of the 800-nm 
beam, which we focus at a few 10 13 W/cm 2 in the jet. 
The probe ellipticity e is controlled by rotating a half 
waveplate in front of a fixed quarter waveplate in order 
to keep the main axis of the polarization ellipse fixed, 
e is defined through F(t) = F cos(u)t)x + ieF sm(ut)y 
and the total intensity |F| 2 is kept constant throughout 
the experiment. We perform this scan for alignment an- 
gles of and 90 degrees. The harmonics are dispersed 
by a 1200 mm^ 1 grating and imaged on microchannel 
plates. The cylindrical geometry of the grating, com- 
bined with an optimized adjustment of the position of 
the gas jet with respect to the laser focus, enables us to 
discriminate between short and long trajectories, spec- 
trally and spatially, in both below- and above-threshold 
regimes [ill, E3] • In the following, we will focus on the 
spectral range in the vicinity of the ionization potentials 
of N2 and O2, which respectively correspond to 10.05 uiq 
and 8.78 wo, where wq is the fundamental laser frequency. 

Figure 1 shows the harmonic signal as a function of 
the laser ellipticity for O2 molecules for alignment an- 
gles of (left) and 90 (right) degrees. Long trajectories 
were measured down to harmonic 11. As we increase the 
laser ellipticity, all the long trajectories' signals (bottom 
panels) decrease significantly. Such a response has been 
explained by the SFA for above-threshold harmonics [l8[ : 
the ellipticity of the driving field induces a lateral shift of 
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FIG. 2: Normalized harmonic yields in rotationally excited 
N2 for harmonics 9-13,19, as a function of the pump-probe 
delay. The solid and dashed lines correspond to the short and 
long trajectories, respectively. 



the free electron wavepacket, and therefore reduces the 
recombination probability. In the vicinity of the ioniza- 
tion threshold, we observe a similar effect. 

As we switch from long to short trajectories, a striking 
difference appears. For an alignment angle of degrees, 
harmonics 9 and 11 increase as a function of the ellipticity 
and reach a maximum at e ~ 0.2 and 0.15, respectively. 
The signal at this point is enhanced by a factor of 9. 8 ±0.7 
for the 9 th harmonic and by 1.5 ± 0.1 for the 11 th har- 
monic. A similar effect is observed for an alignment angle 
of 90 degrees, but only for the 9 th harmonic (2.7 ± 0.2). 
When the intensity / is slightly increased (from ~ 1.6 
to ~ 1.8 x 10 14 W/cm 2 ) this enhancement is observed 
in an additional harmonic - 13 and 11 for and 90 de- 
grees, respectively. Well above the ionization threshold, 
all short trajectories responses exhibit the usual decaying 
rate, regardless of the intensity. 

Repeating the experiment with N2 molecules, no 
anomalous responses are observed. All harmonic orders 
decay smoothly as a function of e, exhibiting similar de- 
cay rates for both short and long trajectories (not shown). 
So, while all long trajectories follow the expected ellip- 
ticity response as predicted by the intuitive three-step 
model, the reactions of near-threshold short trajectories 
are system-dependent and highly selective. The observed 
enhancements at e 7^ cannot be related to the structure 
of the ground molecular orbital, and related interference 
processes at the time of recombination since these 
latter would affect both short and long trajectories. We 
shall return to this point later on, in the context of our 
theoretical analysis. 

An additional difference between short and long tra- 
jectories in N 2 molecules near I p is observed when the 
pump-probe delay is scanned using linearly polarized 
beams. Figure 2 shows the results of our pump-probe 
scans around the half-revival of N2 [16], for both short 
and long trajectories. As observed in previous experi- 
ments, the total harmonic signal peaks at a delay of 4.1 
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ps, when the molecules are aligned along the polarization 
axis of the pump pulse, and has a minimum at 4.3 ps, 
when the molecules are anti-aligned (see Fig. 2 a,c, and 
d for harmonics 9, 13, and 19). Long trajectories were 
measured down to the 9 th harmonic, with the typical re- 
vival structure (dashed lines in Fig. 2). However, the 
11 th harmonic has surprising features. Whereas the long 
trajectory exhibits the typical revival shape, with con- 
trast ™ a l-™™ = 0.10 ± 0.03, the short one (solid line) 

0.35 ±0.04). 
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In order to verify that the inversion effect we observe does 
not depend on propagation, we performed a systematic 
study by varying the gas pressure by one order of mag- 
nitude and found similar results. An intensity scan from 
1 to 2xl0 14 W/cm 2 also shows the robustness of this ef- 
fect. Only below 1.2xl0 14 W/cm 2 does the inversion of 
harmonic 11 disappear. 

Both experiments - the cllipticity and revival scans - 
indicate that threshold harmonics are generated by two 
completely different processes. Whereas the long trajec- 
tories behave as higher harmonic orders and therefore 
seem to be generated according to the standard three- 
step model, it is clear that the short trajectories orig- 
inate from a different mechanism. How can threshold 
harmonics be explained in the framework of classical tra- 
jectories? What breaks the symmetry between the two 
families of electron trajectories? 

In order to answer these fundamental questions, we 
have performed statistical Classical Trajectory Monte 
Carlo (CTMC) calculations 20] that treat the laser V F 
and Coulomb Vq potentials on the same footing. We 
consider the prototypical N 2 system described in terms 
of a one-dimensional scaled hydrogen atom with an ef- 
fective nuclear charge Z e such that Z e — y/2 I p . The 
CTMC procedure employs an N-point discrete repre- 
sentation of the phase-space distribution, g(x,p,t) = 

jj Ylj=i $( x ~ x j{t))${P ~ Pj(t))i m terms of N non- 
interacting trajectories {xj(t),pj(t)}j = i^^N. The initial 
condition, g(x,p,t = 0), consists of the Wigner phase- 
space distribution associated with the fundamental clas- 
sical energy bin, as detailed in [20l ]. 

The CTMC picture of HHG is obtained by consider- 
ing an instantaneous recombination to the ground state 
at time t' whenever a trajectory with energy E(t') — 
p(t') 2 /2 + Vc{t') returns to x(t') = after having been 
ionized at time t" . Such a process is accompanied by 
the emission of a photon with energy E(t') + I p . Fig- 
ure 3 presents the harmonic order of the photon emitted 
upon recombination as a function of time t' of emission. 
Our CTMC results are compared to their SFA counter- 
parts. For above-threshold emission, CTMC and SFA 
results are in good agreement. In contrast to the SFA, 
the CTMC method yields long trajectories for below- 
threshold harmonics; in the (£',H) plane of Fig. 3, 
these trajectories appear as continuations of the above- 
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FIG. 3: Calculated harmonic order as a function of time of 
emission t' , for N2 embedded in a two-cycle field with A = 800 
nm and I = 1.5 10 14 W/cm 2 . The blue (empty) circles 
correspond to SFA calculations, whereas black (dark) and 
green (light) dots, respectively, represent harmonics stemming 
from ionizing and excitation trajectories, in the framework of 
CTMC-Wigner calculations. 
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FIG. 4: Calculated electron position x (a,c), and energy E 
(b,d) as functions of time t, for some short (right column) 
and long (left column) trajectories of Fig. 3. The dotted line 
in (a,c) represents the evolution of the electric field F(t). 



threhold ones below I p . 

To illustrate how these trajectories succeed in return- 
ing to the ion with E(t') < 0, we present in Fig. 4 the 
temporal evolution of the electron position x(t) and of 
the energy E{t) = p(t) 2 /2 + Vc(t), for some typical long 
trajectories. Ionization occurs at t" ~ 0.3 To, where Tq 
is the laser period, soon after the field maximum. At this 
time, the dispersion of the free wavepacket in coordinate 
space, although small, is sufficient to allow the subse- 
quent differentiation of the trajectories in the combined 
laser and Coulomb fields. The electron returns to the ori- 
gin x = at 0.95Tb < t" < 1.3T . Earlier electrons emit 
higher harmonics. At the time of recollision, the laser 
force —F(t) acting on the electron is in the opposite di- 



rection to the Coulomb force [13j. Negative energy E(t') 
is enabled for those later trajectories that recombine at 
t' ~ 1.25T , where p 2 {t')/2 + V C (f) < 0. Our CTMC 
model thus shows clear continuity between above- and 
below-threshold harmonics for the long trajectories. All 
of them obey the standard three-step process and there- 
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fore have similar responses to ellipticity and revival scans. 

When the same analysis is applied to short trajectories, 
no below-threshold harmonics appear (see Fig. 3). We 
thus modify our model by relaxing the condition regard- 
ing ionization: we only require the trajectory to leave the 
fundamental classical energy bin before it returns back 
to x = 0, but we do not require it to be ionized (see 
Fig. 4). In other words, we allow excitation pathways 
to contribute to HHG. The results corresponding to this 
new criterion are included in Fig. 3; below-threshold 
short trajectories emerge, which is consistent with our 
experimental observations. They appear as the contin- 
uation below I p of the above-threshold short trajecto- 
ries. Their temporal evolution is illustrated in Fig. 4(c- 
d). The short trajectories corresponding to u> > I p are 
ionized shortly before recombination, whereas for below- 
threshold HHG, the trajectories recombine before ion- 
ization, with E(t') < 0. Inversely to what occurs for 
long trajectories, earlier electrons emit lower harmonics. 
Clearly, part of the below-threshold HHG signal stems 
from excitation processes, which are classically described 
in terms of short trajectories. Interestingly, a Bohmian 
description of the interaction yields quantum trajectories 
that nicely match, in both position and energy scales, the 
short CTMC trajectories [21|. Those trajectories do not 
involve tunnel-ionization, in agreement with the semi- 
classical calculations of [15j, and therefore are beyond 
the scope of the conventional three-step model. 

We can now reinterpret the experimental observations 
concerning short trajectory below-threshold harmonics in 
light of their excitation origin. With respect to the el- 
lipticity dependence observed in Fig. 1, a description 
of excitation in terms of a non-resonant absorption of H 
photons leads to harmonic yields that smoothly decrease 
as e increases, in agreement with the N2 shapes. This reg- 
ular behavior is modified if intermediate resonances come 
into play during the multiphoton absorption, as shown in 
22] for atomic targets. As e increases, the interplay be- 
tween the decreasing linear and increasing circular field 
components then yields e-responses similar to those of 
77 = 9 — 13 in O2. Obviously, the occurrence of reso- 
nances depends not only on the target considered, but 
also on the alignment angle that determines the allowed 
symmetry of the intermediate excited states. Resonantly 
enhanced excitation may also induce the inversion of re- 
vival structures. Depending on the symmetries of the in- 
termediate states involved in the multiphotonic pathway, 
the harmonic strength can be maximual for a different 
molecular alignment than in the conventional three-step 
process involving tunneling. Finally, it is well known 
in strong-field physics that the ionization threshold of 
a target embedded in an intense laser field is subject 
to an ac-Stark ascending shift proportional to F 2 . In 
this respect, the short trajectories associated with har- 
monic orders close to the unperturbed I p recombine at 
time t', where F(t') ^ 0. Therefore, harmonic orders 



77 > Ip/oJo, which belong to above-threshold HHG at 
low 7, fall down in the below-threshold regime as 7 in- 
creases. Accordingly, we understand the appearance of 
additional anomalous responses to e in O2 as 7 increases, 
and why the inversion of harmonic 11 disappeared in N2 
for 7 < 1.2 x 10 14 W/cm 2 , where 77 = 11 then effectively 
corresponds to above-thrchold HHG. 

To sum up, our joint experimental/theoretical study of 
harmonic generation in aligned N2 and O2 has allowed 
us to observe and explain, unambiguously, two main con- 
tributions to near-threshold HHG. The long electron tra- 
jectories belong to the three-step model, and encode the 
structural information in the same way as in conventional 
high harmonic spectroscopy. The short trajectories stem 
from multiphoton-driven pathways, and can be used to 
reveal complementary information on excited states of 
the molecule. The ellipticity and alignment responses of 
the short trajectory signal will provide information on 
the symmetry and dynamics of the intermediate states 
involved in the multiphoton process. This approach adds 
to the ability of HHG spectroscopy to follow a molecular 
orbital during chemical dynamics 23] . Since this method 
is extremely sensitive to the location of I p , it is possible 
to probe changes in the threshold itself during chemical 
processes. Thus, our study is an integral part of the ongo- 
ing effort to resolve the structure as well as the dynamics 
of molecular orbitals. 
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